Geometrical phase effects in biaxial nanomagnetic particles 
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The oscillation of tunnel splitting are obtained by geometrical analysis of the topological Wess- 
Zumino phase on the basis of tunneling paths in biaxial nanomagnetic particles with magnetic field 
along the hard anisotropy axis. This theory not only just yields the previous quantum interference 
results for the ground state tunneling, but also gives the excited level splittings, all of which agree 
well with the numerical diagonalization. Furthermore, the parity effect in the asymmetric double well 
system which was recently discovered in experiment, can be derived by similar arguments, and is also 
certified by using the complex periodic orbit theory. The possibility of improving the discrepancy 
with the experiment in the periods of oscillations is discussed. 
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I. INTRODUCTION 

For last decade, many physicists have explored the 
nanoscale magnetic particles to study the extrapolation 
of quantum mechanics to macroscopic realm. The mag- 
netizations of these particles show the macroscopic quan- 
tum tunneling (MQT) from a metastable state to a stable 
one, or macroscopic quantum coherence(MQC) between 
two degenerate ground states separated by potential en- 
ergy barrier 0-^]. In the latter case, the geometrical 
phase, known as the Wess-Zumino phase Q|, plays a cru- 
cial role in tunnel splitting. In particular, the biaxial 
nanomagnetic particles present some exotic phenomena 
related to this topological phase. In the absence of exter- 
nal magnetic field, the ground state tunneling vanishes 
only for half-integer spins, which is called as spin par- 
ity effect and related to well-known Kramer degeneracy 
||. When the magnetic field is applied along the hard 
anisotropy axis, however, the tunnel splitting oscillates 
with magnetic field regardless of spin values ||. These 
spin parity effect and oscillations of tunnel splittings are 
interpreted as the interferences of geometrical phases be- 
tween two opposite tunneling paths in the instantonic 
approach [|| [)| . 

Most theoretical treatments so far about the oscilla- 
tions of level splittings have been restricted to the in- 
stanton method for the mapped effective one-dimensional 
potential by performing the cos integration in spin 
coherent state representation |?],|8| or by direct particle 
mapping method In these approaches, although the 
ground state tunneling can be obtained well, it is not 
possible to treat the tunneling between excited energy 
levels exactly. In fact, the result of numerical diago- 
nalization for excited states is that the number of os- 
cillations decreases one by one as the quantum number 



increases. In addition, due to the mapping onto one- 
dimension, the dominant tunneling path between two 
wells in (9, (j>) space is not clear. This makes it diffi- 
cult to perform the direct geometrical analysis of Wess- 
Zumino phase. Recently, Garg jl(J derived the splitting 
oscillations in high energy levels using the discrete WKB 
approximation. Another semiclassical approach is the 

the extension of 



complex periodic orbit theory |11 
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the trace formula in periodic orbit theory [|1 3| — 11 5f| to tun- 
neling systems. We obtained successfully the energy lev- 
els for not only the ground state but also all the higher 
excited states by applying this method to spin system 
for the first time, and showed the complete agreement 
with numerical diagonalization [fl6|| . Since this theory 
also uses the effective mapping potential obtained by cos 
^-integration of original Hamiltonian, the dominant tun- 
neling path in phase space remains to be unspecified yet. 
The level splitting oscillations were confirmed by experi- 
ment very recently in the Fes molecular clusters with spin 
S = 10 |0. They showed that, although the oscillations 
occurrs as predicted by theory, in quantitative aspect the 
oscillation periods is about 1.6 times larger than previous 
theoretical results, and this discrepancy can be resolved 
introducing the fourth-order terms in spin Hamiltonian. 

In this paper, we find the dominant tunneling paths in 
phase space, and analyze geometrically the Wess-Zumino 
phase without any mapping in order to obtain the ground 
and excited level splitting oscillations with magnetic field 
along the hard axis. Furthermore, the experimentally 
discovered parity effect jl7[| which is very similar to spin 
parity effect in the asymmetric system is derived theo- 
retically using this geometrical analysis and complex pe- 
riodic orbit theory. We also discuss the possibility of 
improving the discrepancy with experiments by adding 
the fourth-order terms in spin Hamiltonian. In section 
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II, after rinding the dominant tunneling path, we obtain 
the splitting oscillations for the tunneling between both 
the ground and the excited states in symmetric system by 
investigating the relation between the geometrical phase 
and the tunneling path in phase space without mapping 
onto a particle. In Section III, the parity effect occuring 
in the asymmetric system is derived using the arguemcnt 
of Section II, and furthermore is also certified from the 
complex periodic orbit theory. We discuss the possibility 
of improving the quantitative discrepancy with experi- 
ment and give conclusions in Section IV. 



II. LEVEL SPLITTING OSCILLATIONS IN THE 
SYMMETRIC SYSTEM 

In this section, we calculate the both ground and ex- 
cited level splitting oscillations by analyzing the topo- 
logical Wess-Zumino phase geometrically in phase space. 
The tunneling rate in spin system is 



T = J V[cos6}V[(t)}exp(-SE) 



(2.1) 



in the path integral formalism, where T> means the in- 
tegration over all paths and Se is the Euclidean action 
which is given as 



Se — iSwz 



dTH{9, . 



(2.2) 



Swz = S / dr(j} [1 - cos 6{t)] 
#{S-*#(t)]} 



Here, r is the imaginary time, H.(0, <fi) is the spin Hamil- 
tonian in the (0, </>) representation, and 



(2.3) 



is the Wess-Zumino action with total spin number S. 
Here p = ScosO and corresponds to the conjugate mo- 
mentum to coordinate 0, since 4> and p satisfy the Hamil- 
ton's equations of motion in the spin system. When we 
consider the tunneling from <j) — to 7r, Eq. (2.3) means 
the area surrounded by p = S, <j) = 0, 7r lines and the 
tunneling path p(4>) in (</>,p) phase space. [See the Fig. |l| 
(a).] Since the energy splitting due to tunneling, AE, is 
determined by Wess-Zumino action (the area mentioned 
above), i.e., 



AE ~ cos(SWz), 



(2.4) 



whenever this area becomes (n + 1/2)tt with integer n 
the tunnel splitting vanishes. In order to perform the ge- 
ometrical evaluation of Wess-Zumino phase without map- 
ping, it is essential to find the dominant tunneling path 
p(<p) in the phase space. It will be shown below that 



in the case of fcth state, the Wess-Zumino phase contin- 
uously decreases from (S — k)n to zero as the external 
field h increases, so that the moment of area being the 
half-integral multiple of n, at which the tunnel splitting 
quenches, occurs (S — k) times. 

Let us consider the Hamiltonian describing the biaxial 
spin system which is given by 



H = -DSl + E(S 2 - S 2 ) + c?mbH • S, 



(2.5) 



where D and E are longitudinal and transverse 
anisotropy constants, respectively, g is the gyromagnetic 
ratio and /j,b is the Bohr magneton. This Hamiltonian 
represents that the system has an easy axis in the re- 
direction, and hard axis in the z-direction. When the 
magnetic field H is applied along the z-direction, the sys- 
tem is symmetric and the reduced Hamiltonian is given 
in the spin coherent state representation by 

H' = [2(1 - A) + (2A - 1) cos 2 <f\p 2 - 2Shp 

— (1 - X)S 2 — (2A — 1)5 2 cos 2 <j), (2.6) 

where H' = H/(D+E), A = D/(D+E) and = H/H c = 
gp-BH /2(D + E)S with coersive field H c . In (p, <fi) phase 
space, li! has two degenerate minima when p = Sh and 
4> = 0, 7r. At 4> = \i however, p — Sh/h\ with hi = 
2(1 — A) gives the saddle point for < h < hi, local 
minimum for hi < h < hi with hi = y2(l — A) and 
global minimum for /12 < h < 1, respectively, if p can 
be also defined beyond the range — S < p < S. These 
pictures of various external fields are shown in Fig. [l] as 
the energy contour plots. 

In order to get an information about the dominant 
tunneling path, we revisit, for a moment, the mapping 
onto a one-dimensional particle problem by integration 
over cos 8 in Eq. (2.1). The resulting Euclidean action is 
given by 



S E (4>) = iA{4>) + SV2X-1 / dr 



\m(^ 2 



V{<t>) , 
(2.7) 



where the dot means the Euclidean time derivative and 
the imaginary part of Euclidean action which is respon- 
sible for the phase effect is 



A{4>) = s 



1 - 



1 - (2A - 1) sin 2 , 
the ^-dependent effective mass is 

MU) = 5— 

y J 1 - (2A - 1) sin 2 (/} 

and the effective potential is 



y(0) = isin 2 



1 



1 - (2A - 1) sin 2 



(2.8) 



(2.9) 



(2.10) 
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For the ground state tunneling, if we compare Eq. ( |2.3| ) 
and Eq. (2.8), the tunneling path p((f>) is found to be the 
second 

It is emphasized that the tunneling path p((f>) in 
the phase space {4>,p) is just the path satisfying 
dH'((j),p)/dp = 0. This means that the tunneling path 
follows the minimum energy positions at all <j) values. 
These tunneling paths for vaious fields are shown in Fig. 
1 as thick solid lines. This simple condition for domi- 
nant tunneling path in the phase space (</>, p) works well 
only when the Gaussian integration is possible as in the 
present case. 

Now, we calculate the Wess-Zumino action (the area) 
mentioned above as h increases from to 1. When h = 0, 
the phase is just Sir (whole area of phase space consid- 
ered) because the dominant tunneling path is p(4>) = 0. 
It is important to note that, although the quantum me- 
chanical ground state energy is not the well minimum, 
the path starts and ends at well minima, because the 
Euclidean action with the energy of the minimum gives 
the information about the ground state splitting in semi- 
classical theory. When < h < hi [Fig. |l| (a)], the path 
is confined within a range < p < S. Therefore, we 
can easily calculate the area enclosed by the path and 
p = S line. Until h becomes hi, where the north-pole be- 
comes the saddle point of energy, the area continuously 
decreases from Sir. For h± < h < h 2 [Fig. |l| (b)], how- 
ever, the local minimum exists outside the p = S line. 
In this case the tunneling path runs beyond the p = S 
line and thus, the Wess-Zumino action becomes the area 
of (region B + region C - region A). As h increases the 
action continues to decrease and becomes zero at h = hi . 
As a result, within a range < h < h 2 , the quench- 
ing take places S times for the ground state by Eq. ( p.4[^ 
The range hi < h < 1 is more or less subtle [Fig! [y 
(c)]. The local minimum turns into the global one, i.e., 
the point {Sh/h\, 7r/2) becomes lower than the well min- 
imum. In this case, the path giving the Wess-Zumino 
action should be determined by somewhat different way, 
because there exists the contour of Eq (well minimum 
energy). If the path obtained by the previous method 
encounters the contour of Eq, it should follow the con- 
tour in a real time due to energy conservation. There- 
fore, the action is equivalent to the area of [region B + 
region C - region A ] of Fig. 1(c) which always van- 
ishes in the range of h 2 < h < 1. This vanishing of the 
action is closely related to the cancellation between res- 
onant tunneling and quenching phenomena fl6| ] and also 
between the Euclidean action of real time motion and 
Wess-Zumino action M. 

The extention to the excited states is straightforward. 
The fcth energy level of each well can be calculated from 
the well-known EBK (Einstein-Brillouin-Keller) quanti- 
zation rule [181, 



p(cj))d(j) = 2?r [k 



k = 0,1,2, 



(2.11) 



In this quantization, whenever the action for closed orbit 
in one well differs by 2-7r, the energy level is given. It will 
be shown that, by this fact, the number of splitting os- 
cillations decreases one by one with increasing quantum 
number. When h = 0, there exist the contours of E^ cen- 
tered at p = and (j> = 0, ir whose areas are 2kw for kth 
excited state (see also Fig. ^). If we use the analogy for 
the case hi < h < 1 of the ground state tunneling, the 
area of kir is excluded, due to the energy conservation, 
for the fcth excited state (kir/2 for each well) in the phase 
space considered, i.e., (S—k)n. As h increases the contin- 
uous reduction to zero is same as the ground state case. 
This leads to the (S — k) oscillations for fcth energy level. 
The calculated periods of oscillations for both ground and 
excited state levels completely agree with the numerical 
diagonalizations. Comparing with the previous semiclas- 
sical theories , in this theory the dominant tunneling 
path in {<f>, p) phase space can be found to make possible 
to perform the geometrical analysis of the effect of topo- 
logical phase, and in addition to the gound state case the 
structures of splitting oscillations for the excited states 
can be obtained easily on the basis of EBK quantization 
scheme. 



III. PARITY EFFECT IN THE ASYMMETRIC 
SYSTEM 

Unlike the symmetric case in Sec. II, in this section 
we treat the asymmetric system having the longitudinal 
field as well as the transverse field. In such system, we 
derive the parity effect which was discovered in recent 
experiment [fL7| , by using the simple geometrical argu- 
ment developed in section II, and also lead to the same 
result in the semiclassical complex periodic orbit theory. 
This parity effect is very similar to the spin parity ef- 
fect proposed by Loss et. al. and von Delft and Henley 
Q. This is observed when one measures the transtion 
between the ground state in upper well and the fcth ex- 
cited state in lower well, by applying the external field h' 
along the easy axis (longitudinal field) which makes the 
system asymmetric in addtion to h in the hard direction 
(tansverse field). For the case of the tunneling from the 
ground state of the upper well to even fcth state of lower 
well, the phase of the oscillation is same as that of the 
symmetric case, while for the case of tunneling to odd 
A;th state of the lower well the phase is shifted by it/2. 

This interesting observation can naturally arises in our 
above analysis. As shown in Fig. 2, the energy structure 
of spin system is asymmetric due to h' . At certain val- 
ues of hi the ground state energy level E l Q in upper (left) 
well exactly coincides with the fcth excited level E r k in 
lower (right) well, so that one can consider the tunnel 
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splittings. The tunneling path starts from the minimum 
at upper well, and when it exits the barrier it follows 
the contour of fcth excited state in the lower well due 
to energy conservation. Therefore, by EBK quantization 
relation and analogous arguments with the excited states 
in Section II, the area kir/2 in the {<f>,p) phase space is 
excluded when h — case. As the transverse field h in- 
creases, total area surrounded by this path continuously 
decreases from (Sir — kir/2) to zero in the same man- 
ner with the symmetric case except the kir/2 shift. This 
leads to the parity effect in the asymmetric system. 

Next, let us consider this effect in complex periodic 
orbit theory. In quantum mechanics the energy spectrum 
of the system can be derived from the singularities of the 
trace of the energy-dependent Green's function which is 
given by 

\ - d„ 



g{E) =Tr (E - H 



E — E n 



(3.1) 



where Tr means the trace and d n is the degeneracy factor. 
The semiclassical approximation of Eq. ( |3.l| ) was given 
by Gutzwiller |l3( ] in the form 



g sc (E) 



E 
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(3.2) 



where Tj, Sj and fj,j are the period, the action and the 
Maslov index of the jth primitive classical periodic or- 
bit, and index r corresponds to the repetition of prim- 
itive periodic orbit. In evaluating this sum, we restrict 
on the one-dimensional effective particle problem. The 
system with b! = which are given as Eq. (2.7) - Eq. 
(2.10) is recently treated in detail within the complex pe- 
riodic orbit theory in Ref. |L(|. When the longitudinal 
magnetic field h! is applied, the effective potential of Eq. 
(2.10) becomes asymmetric as shown in Fig. 3. This 
figure presents the case when the ground state energy of 
upper well coincides with the fcth excited level of lower 
well. The classical orbits in two wells and tunneling or- 
bits within the barrier are represented as solid and dotted 
lines, respectively. The directions denoted by arrows are 
given by the rules of Ref. |IlJ . 

Now we should consider all possible periodic orbits in 
order to perform the summation in Eq. (3.2). Let us 



take the classical segments in the respective wells as the 
half of whole orbit having the actions W\ = S\/2 and 
Wi = 6*2/2, and similarly for the tunneling segments 
within the barrier having an action = S c /2. Let us 
consider the case where the ground state energy in the 
upper well is equal to the fcth excited state in the lower 
well. Note, then, that the relation between two actions 
W\ and W2 is 



W 2 = W 1 + kir, k = 0,1,2, ■ 



(3.3) 



due to EBK quantization rule. The possible periodic 
orbits can be parametrized as follows: i) the paths that 



start and end at 4>a or 4>A' without complete rotations, 
ii) the paths that start at ^>a(^>A') and end at <j>A'(<l>A), 
i.e., having complete rotations (See Fig. 3. (f>A and cj>A' 
are the equivalent position in <j) coordinate). The Wess- 
Zumino phase contributes only to the paths belonging to 
category ii) as 



2vrS* ( 1 - — 
h 2 



(3.4) 



For the paths that start at 4>A and end at cf>A', i-e., pro- 
ceed from left to right, we take the Wess-Zumino phase 
as positive, whereas for the paths from right ((/>a>) to 
left (4>a) we take as negative. The Maslov index is eas- 
ily evaluated following the method for double well po- 
tential in Ref. jjlljj , i.e., each classical segment Wi, W% 
accompanies the additional phase +7r/2, while the tun- 
neling segment —tt/2. The all possible primitive peri- 
odic orbits are drawn graphically in Fig. |^, where the 
solid line corresponds to exp[i(Wi-|-7r/2)], thick solid line 
to exp[i(W2 + 7r/2)] and dotted line to exp[— 6 — in/2]. 
Therefore, the semiclassical trace of Green's function can 
be expressed as the summation over all repetitions of 
these orbits, and its poles are found to be when 



I _ e i(2Wi-Hr) 
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l_ e i(2W 2 + x) +e -2e 



e i(W!+W 2 +^) ( e ict 



i(2Wi- 



i(2W 2 +-ir) 



= 



(3.5) 



If we neglect, for the moment, the tunneling probability 
, then the poles give the EBK energies E n for upper 
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well where 



S 1 (E n )=2W 1 (E n )=2w n 



(3.6) 



with n = for the ground state. Therefore, we can ex- 
pand the actions around the EBK energy like 



2W\ (E) = ir + Ti(E — Eq) + ■ ■ ■ , 
2W 2 (E) =Tr + 2kir + T 2 (E - E ) 



(3.7) 



near the ground state of upper well, where T\ and T 2 
are given by dS\/dE and dS 2 /dE, respectively. Here, we 
used the fact that the difference of two actions at EBK 
energy is 2irk for the transition between ground state in 
uppe r we ll and fcth excited state in lower well. Expanding 
Eq. (3.5) up to the first-order in e _2e , 



T{F 2 {E - Eo) 2 - e~ 2e [(-l) fe (e M + e -' la ) + 2] = 



(3.8) 



The final energy splitting is obtained as 
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E = E ± 

E = Eq zb 



2 



cos — , even k, 



-==e- e sin-, oddfc, (3.9) 



which is just the parity effect in the asymmetric system 
of Ref. @. 



IV. DISCUSSIONS AND CONCLUSIONS 

Until now, we have analyzed the Wess-Zumino phase 
of biaxial spin system in (p, (j)) phase space, and obtained 
not only the same results with those in the mapping for- 
malism for ground state, but also the excited energy level 
splitting oscillations. Furthermore, the parity effect in 
the asymmetric system shown in recent experiment |l7]] 
is derived by both the geometrical analysis and the com- 
plex periodic orbit theory. However, the quantitative dis- 
crepancy with experiment in the period of oscillation ex- 
ists, i.e., about 1.6 times larger oscillation period than 
expected in the theory. In order to resolve these dis- 
crepancies the authors of Ref. introduced the fourth- 
order term C(S\ + 5*1) in spin Hamiltonian (Eq. ( |2.0| )), 
where C is the adjustable parameter which is shown to 
be —2.9 x 10~ 5 K in Kelvin unit through numerical di- 
agonalization, when anisotropy constants are given by 
D = 0.292K and E = 0.046K. Then the fourth-order 
Hamiltonian is expressed in the ^-representation as 



Hi=H' + 2C(S* + S*- 6S 2 z Sl) 



(4.1) 



within the semiclassical approximation. In fact, if we 
consider the classical commutation relation satisfied by 
angular momenta, the cubic and quadratic terms must 
be added in this Hamiltonian. But we have ignored them 
since their contributions are negligible. It is not possible 
to get an anaytical effective potential through integra- 
tion, since Eq. ( [4.l|) is beyond the Gaussian approxima- 
tion. However, the qualitative analysis of this Hamilto- 
nian can be performed within our theory on the basis of 
the energy structure of "Hi and the dominant tunneling 
paths in ((f), p) phase space. 

The structure of energy barrier for Tlx is different from 
that for TC as shown in Fig. [| Since C is negative, the 
effect of SlSy in Eq. (1.1) is to raise the barrier around 



7r/2 and the mid- value of < p < S, whereas S^(Sy) 
to lower the barrier at z(y)-axis. Therefore, as the field 
h increases, the saddle point moves more slowly to the 
north-pole, compared with the C — case. If we use 
the same analogy about finding the dominant tunneling 
path, this means that the slower decrease of the area sur- 
rounded by tunneling path, and thus the larger oscillation 
period. Although it is possible to get a qualitative correc- 
tion to C = case in the desirable direction, in order to 
be consistent with experiment, the adjustable parameter 
C must be about 3 times greater than the present value. 



The quantitative disagreement implies that the tunneling 
paths which follow all energy minima are no longer the 
dominant paths in the quartic case. In order to treat the 
quartic case exactly, new theoretical approach in order 
to obtain the dominant tunneling path in phase space is 
needed. 

In conclusion, we derive both the ground and excited 
tunnel splitting oscillations in the biaxial nanomagnetic 
particle with the magnetic field along the hard anisotropy 
axis by finding the dominant tunneling path and geomet- 
rically analyzing the topological Wess-Zumino phase in 
the phase space. All the results of analysis are in agree- 
ment with the numerical diagonalization. Furthermore, 
the interesting parity effect in the asymmetric case is 
naturally clear in this analysis, and also certified within 
the complex periodic orbit theory. We also discussed the 
possibility of improving the discrepancies with experi- 
ment in the period, by introducing the quartic terms in 
spin variables into Hamiltonian. In order to resolve the 
discrepancy quantitatively, the new approach to find the 
dominant tunneling path in quartic case is required. 
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FIG. 1. The contour plots of Hamiltonin Eq. ( 2.6) and the 
tunneling paths (thick solid lines) which follow all the minima 
with p for various external magnetic fields h. Here, A = 0.861 
of Feg is used. In this case, hi — 0.278 and h,2 = 0.528. 
The dashed lines mean p = S lines, a) h — 0.2. The area 
corresponding to the geometrical phase is the region between 
the tunneling path and p — S line, (b) h = 0.5. The phase 
is evaluated as the area of (rigion B + region C - region A), 
(c) h = 0.7. When the path encounters the contour of Eq, it 
follows the contour in real time (thick solid line with arrow). 
The phase is also the area of (region B + region C - region 
A) but yields zero. For the excited states, see the text. 



FIG. 2. The typical asymmetric contour plot of the Hamil- 
tonian with both longitudinal (h!) and transverse (h) mag- 
netic fields. The dominant path starts at the left well mini- 
mum E , and when it encounters the contour of fcth excited 
state in right well i?£(= Eg), it follows the contour (thick 
solid line with arrows). 



FIG. 3. The asymmetric effective potential when the lon- 
gitudinal field is applied. The actions are denoted by Si, S2 
and S c . Actually, 4>a and <j>A' are same positions, but we use 
the different symbols in order to distinguish the direction. 



FIG. 4. The diagram representation of all possible peri- 
odic orbits in the asymmetric potential. The segments repre- 
sent e^i+Vz) ( soM line), e i(w 2 +^/2) ( - tMck goM ^ and 
e -e-i*/2 ( dotted ]ine ) 



FIG. 5. The energy barrier of Hamiltonians Tii(p, 4> = t/2) 
for C = K (solid line), and C = -2.9 x 10~ 5 K (dashed line) 
at h = 0.15. As h increases, the barrier minimum (saddle 
point) moves to the right. But the movement of that for 
C — K case is faster than for C = 2.9 x 10~ 5 K case. If we 
can think these minima as the passage that is taken by the 
dominant tunneling path, the area reduction in phase space 
is smaller for C — 2.9 x 10~ 5 K case. 
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